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Explicit expressions for the concurrence of all positive and trace-preserving ( "stochastic" ) 1-qubit 
maps are presented. By a new method we find the relevant convex roof pattern. We conclude that 
two component optimal decompositions always exist. 

Our results can be transferred to 2 x n-quantum systems providing the concurrence for all rank 
two density operators as well as a lower bound for their entanglement of formation. 
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INTRODUCTION 

In quantum physics a system in a pure state 7r = \ip) (ip\ 
may have subsystems in states which are not pure but 
mixed. These mixed substates are typically correlated in 
a non-local and non-classical way. The use of this phe- 
nomenon of entanglement as a resource for communica- 
tion and computation is a main feature of quantum infor- 
mation theory [l| . This makes the search for a quantita- 
tive understandingand characterization of entanglement 
a central issue [1,13 ■ Entanglement measures ought to 
describe single- use or asymptotic capabilities of quantum 
systems and channels just as the von Neumann entropy 
S(p) = — Trplogp is an asymptotic measure for infor- 
mation content. They are, similar to entropy, non-linear 
and unitarily invariant functions on the space of states. 

Bennett et al. 0] introduced the entanglement of for- 
mation Ep(p) as the asymptotic number of ebits (max- 
imally entangled qubit pairs) needed to prepare the en- 
tangled bipartite state p by local operations and classical 
communication (LOCC) and showed that 



E F {p) = min^p 3 S (Tr B (^j)) 



(1) 



where Tr^ is the partial trace over one of the two subsys- 
tems and the minimium is taken over all possible convex 
{J2Pj — ^-iPj > 0) decompositions of the state p into 
pure states 



nj pure. 



(2) 



Closed formulas for the entanglement of formation, i.e., 
analytic solutions to the global optimization problem ([1]) 
are only known for certain classes of highly symmetric 
states [a, 0] and for the case of a pair of qubits (2x2 



system). In the latter case, the analytic formula for the 
entanglement of formation was obtained first for special 
states |4|, |7| and later proved for all states of a qubit pair 
Q. It expresses Ep(p) in terms of another entanglement 
measure C(p) which was named concurrence in [7[. The 
concurrence appeared to be an interesting quantity in 
itself @. 



ToL 11, 13], have obtained 



Many authors, e.g. 
bounds for the concurrence of larger bipartite systems. 

In the present paper we obtain analytic expressions 
for the concurrence for general stochastic 1-qubit maps 
and therefore for general 2 x n bipartite systems provided 
the input state p has rank two. For this we employ the 
convex roof construction 13|, LL4] as a way to study global 
optimization problems of the type ([I]) . Our main results 
are given by Theorems 2 and 3. 

Let $ be a positive and trace-preserving (i. e., stochas- 
tic) map from a general quantum system into a 1-qubit- 
system. This setup includes as special case the partial 
trace Tr^ which maps states of a bipartite 2 x n sys- 
tem to states of the subsystem. For pure input states 
tt = \ip)(ip\ the concurrence is defined as 



C$(7r) = 2v/det$(vr) 
and for a general mixed input state p one defines 



(3) 



(4) 
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where the minimum is again taken over all possible con- 
vex decompositions into pure states. Let us consider the 
case where p has rank 2 and is therefore supported by a 2- 
dimensional input subspace. Then we have to consider in 
^ only pure states supported in the same 2-dimensional 
supporting input space. By unitary equivalence we are 
allowed to identify input and output subspaces. Hence, 
calculating the concurrence of a rank two density opera- 
tor p = j—i Pij \vi) (vj | of a 2 x n system is equivalent to 
computing the concurrence of a certain 1-qubit stochas- 
tic map. This map is completely positive and explicitely 
given by $(p) = Pij D ij witn Aj = Tr B \vi)(vj\. 
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However, our construction of the concurrence works 
for all stochastic 1-qubit maps, not only for completely 
positive ones. It is therefore suggestive, but not the topic 
of the present paper, to ask for applications to the en- 
tanglement witness problem [l5j |. 

In some cases the convex roof for the concurrence ap- 
pears to be a flat convex roof. In these cases optimal 
decompositions for the concurrence also provide optimal 
decompositions for the entanglement of formation and 
therefore Ep(p) can be expressed as a function of the con- 
currence C(p), exactly as in the case of a pair of qubits 
8]. If the roof of the concurrence is not flat, our results 
for the concurrence provide a lower bound for the entan- 
glement of formation. 

We illustrate our procedure by explicit formulas for 
the cases of bistochastic and of axial symmetric 1-qubit 
maps. In both cases the result is of a surprising trans- 
parency. 



THE CONVEX ROOF CONSTRUCTION 

In the following, all linear combinations are under- 
stood as convex combinations, i.e., the {pj} always sat- 
isfy ^2pj = 1 and pj > 0. Solutions to the optimization 
problem eq. (|4l2p can be characterized as so-called convex 
roofs: Let f2 denote the convex set of density operators p 
and let g(jr) be a continuous real- valued function on the 
set of pure states. 



that 



Theorem 1 (see 13, HI)- There exists exactly one func- 
tion G(p) on f2 which can be characterized uniquely by 
each one of the following two properties: 

1. G(p) is the solution of the optimization problem 



G{p)= min VpjstTTj). 

P=2_ J Pi "j — 



(5) 



2. G(p) is convex [16( and a roof i.e., for every p 6 
exists an extremal decomposition p = ^Pjftj such 



G(p) = ^Pjgfa) 



(6) 



Furthermore, given p, the function G is linear on the 
convex hull of all those pure states which appear in the 
decomposition ^ of p. Therefore, G provides a foliation 
of fl into leaves such that a) each leaf is the convex hull 
of some pure states and b) G is linear on each leaf. If G 
is not only linear but even constant on each leaf, it is a 
flat roof. 

STOCHASTIC 1-QUBIT MAPS 



The space Ai 2 of hermitian 2x2 matrices p 
is isomorphic to Minkowski space R 1,3 via 



(x ,x) 



1 



{x Q I 



x ■ a) 



Cll ) 



(7) 



1 I x Q + x 3 x% + ix 2 

2 \xi - ix 2 x - x 3 



We have det p = j(xq — x\ — x\ — x\) = jX-x where the 
dot between 4- vectors denotes the Minkowski space inner 
product and Trp = xq. Therefore the cone of positive 
matrices is just the forward light cone and the state space 

of a qubit, the Bloch ball, is the intersection of this 
cone with the hypersurface V defined by xq = 1. In 
this picture mixed states correspond to time-like vectors 
and pure states to light-like vectors, both normalized to 
x = 1. 

A trace-preserving positive linear map $ : A4 2 — > M 2 
can be parameterized as [17| 



*!>(/;): <I> ( ^(x I + x-a] 



1 



XqI + (xot + Ax) ■ <t) 



(8) 



where A is a 3x3 matrix and t a 3- vector. 

We consider the quadratic form q on M 2 defined by 



g*(x) = 4(det $(p) - w det p) = $(x) • $(x) - w x • x = 



(9) 



where w is some real parameter. For pure states, i.e., 
on the boundary of the Bloch ball where x • x = 0, 
the form q(x) equals the square of the concurrence C = 
2 v / det¥(p). 

Furthermore, we denote by Q the linear map Q : xi 1— > 
^2 QijXj corresponding to the quadratic form q via polar- 



ization: 

~$ fl-\t\ 2 -w -tA \ , n . 

Q w = Q -wr }ij = { _ (fA)T wI _ a t a ) (10) 

The central result of this paper are the following two 
statements: 
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Theorem 2. Let the quadratic form q and therefore the 
matrix Q be positive semidefinite and degenerate, i.e., 
Q > and dim Ker Q > 0. If KerQ contains a non-zero 
vector n which is space-like or light-like, n • n < 0, then 
g 1 / 2 is a convex roof. Furthermore, this roof is flat if such 
an n exists with no = 0. 

Theorem 3. For every positive trace-preserving map <& 
exists a unique value Wq for the parameter w such that 
the conditions of Theorem 2 are fulfilled. Therefore the 
concurrence of an arbitrary stochastic 1-qubit map <& is 

given by C$(p) - 



Let us sketch the proof of Theorems [2] and [3] The 
square root Jq of a positive semidefinite form q on a 
linear space provides a seminorm on this space and is 
therefore convex. According to Theorem Q] we need to 
show that it is also a roof, i.e., there is a foliation of the 
space into leaves such that q 1 ^ 2 is linear on each leaf. Let 
n = (no,n) be a non-zero vector in KerQ. Then for all 
vectors m we have 

q(m + n) = (m + n)Q(m + n) = mQm = q(m). (11) 

Let us start with the case where n can be chosen to have 
no = 0. Then n gives a direction in V along which q is 
constant. Therefore, y/q is a flat convex roof. 





















sr > 



M 2 



FIG. 1: The embedding of the Bloch ball into M2 and its 
foliation by a flat convex roof. 

Let us now consider the case where Ker Q does not con- 
tain a vector n with no — 0. Then we have dim Ker Q = 1 
and this line intersects V in one point which we call n. 
Every other point m in V can be connected to the point 
n by a line lying in V. 
half- line K + 3 sm 4 



Then q 1 / 2 is linear along the 



(1 — s)n since 

q (sm + (1 — s)n) = (sm + (1 — s)n)Q(sm + (1 
= s 2 q(m) 



s)n) 
(12) 



This concludes the proof of Theorem [2l 

For the proof of Theorem [3] we note that the space V 
of stochastic maps is itself a convex space. It can be pa- 
rameterized as follows 13: Let i be a unit 3- vector and 



KerQ 




FIG. 2: The foliation of the Bloch ball in the case no ^ 0. 

a, /?, ui, W2, W3 be parameters taking values between zero 
and one: < a < 1; 0</3< 1; < < uj 2 < lo 3 = 1. 

With the abbreviation v = \Jy^h=i we can repre- 

sent stochastic maps © up to orthogonal transforma- 
tions by t — (ti,t2,t3), A = diag(Ai, A2, A3) where 



U 
A,: 



aw,- 



(13) 
(14) 



Furthermore, the boundary dV is given by j3 = 1. In 
this case, the unit vector £ represents the touching point 
(or one of the touching points in more degenerate cases) 
between the unit sphere and its image. Let <!> £ dV, so 
(3=1. Then it is easy to check that w = av 2 makes Q 
positive semidefinite since it permits a Cholesky decom- 
position Q = RR T into a triangular matrix with a zero 
on the diagonal: 



R = 



\_ 

v[i\ 

vpii 

V-^36"3 v[i 3 ) 



(15) 



where fii = y/ct(l — otw 2 ). Furthermore, n = (1, ^Ci^i) 
is a lightlike vector in KerQ. 

In the general case < 1 we have 



P 2 Q 



w0- 



+ (l-/3 2 ) 



^10 0^ 



\o 0/ 



(16) 



Therefore, Q 



as a sum of two positive semidefinite 
terms is either positive semidefinite or positive definite. 
In the first case we are done with wo = af3 2 v 2 . In the 
other case we must adjust wq ■ It is clear that Q is not pos- 
itive for w — > ±00. Therefore due to continuity, we can 
make Q positive semidefinite and degenerate by increas- 
ing or decreasing w. Let wi < u>2 be the points of de- 
generation and ni,ri2 corresponding vectors in KerQ,^. 
Then (eq. ([ID))) niQ ni = win 2 and n 2 Qon 2 = w 2 n|. 
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Furthermore, no nonzero vector can be both in Ker Q Wl 
and KerQ W2 . So, niQ ni > I02ii| and n 2 Qo n 2 > winij, 
providing (w 2 — Wi)n 2 < and (u> 2 — Wi)n| > 0. There- 
fore, increasing w will make Ker Q time-like and decreas- 
ing w will make it space-like. This proofs the claim of 
Theorem^ existence of a suitable u> . Uniqueness can be 
shown easily. It also follows indirectly from the unique- 
ness of the convex roof extension, Theorem [TJ More de- 
tails can be found in 11911. 



EXPLICIT EXAMPLES 

Let us demonstrate our construction on some exam- 
ples. 

Bistochastic maps or unital channels 

Unital 1-qubit channels are quite trivial. We have t = 
0, A = diag(Ai, A2, A3) and therefore w — max(Af , A|, A3) 
fulfills the conditions of Theorem [2] and provides the roof 



In the case 1 > 1 we have KerQ = Span{e x , e y } and 
the resulting roof is flat. In the other case we have a 
one-dimensional Ker Q generated by n — (1,0, 0, Zq) with 

2q = VTO^+VgE S and a non-flat roof. 
VtI 1 — i')-v a ( 1_Q ) 



CONCLUSION 

We calculated the concurrence C$ of all trace- 
preserving positive 1-qubit maps and therefore for gen- 
eral 2 x n bipartite systems with rank-2 input states. 

The concurrence is real linear on each member of a 
unique bundle of straight lines crossing the Bloch ball. 
The bundle consists either of parallel lines or the lines 
meet at a pure state, or they meet at a point outside 
the Bloch ball. Furthermore, C$ turns out to be the 
restriction of a Hilbert semi-norm to the state space. 

More details and applications, including the entangle- 
ment of formation in 2 x n systems and the Holevo ca- 
pacity [2(| of channels will be given in [u| • 



C{ P )=q X '\p) 



\ 



(l-w)xZ + j2(w->Z)x* (17) [1 



which is flat in one direction since one of the terms in the 
sum vanishes. 

Nevertheless, this case includes channels of all Kraus 
lengths between 1 and 4. 



Axial symmetric channels 

Every channel commuting with rotations about the x%- 
axis is of the form 



Hp) 



ax 00 + (1 - j)xi 



0x 
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f3x 01 

jxn + (1 - a)x 00 



(18) 



This corresponds to A = diag(/3, (3, a + 7 — 1) and 
t = (0,0, a — 7). This family includes many standard 
channels, e.g., 

• the amplitude-damping channel (length 2, non- 
unital) for 7=1, (3 2 = a; 

• the phase-damping channel (length 2, unital) for 
a = 7 = 1 and 

• the depolarizing channel (length 4, unital) for a = 
7, P = 2a-1. 

Here we find that ql/ 2 is a convex roof for w = 
max(/3 2 , 1 ) with 

01 = 1 + 2a7 - a - 7 - 2y/a(l - a) 7 (l - 7). (19) 
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